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ABSTRACT 


In this paper, we study biharmonic isometric immersions of a surface into and 
biharmonic Riemannian submersions from the product space M? x IR. We give 
a classification of proper biharmonic isometric immersions of a surface with 
constant mean curvature into M? xR. More precisely, we prove that a surface 
with constant mean curvature H in M? x R is proper biharmonic if and only 
ifit is a part of S! (zm) x R in S? (at) x R. We also obtain a complete 
classification of proper biharmonic Hopf cylinders in M? x R. On the other 
hand, we give a classification of biharmonic (including harmonic) Riemannian 
submersions 7 : M? x R — (N?,h) from the product space, and construct 


many family of proper biharmonic Riemannian submersions M? x R > R?. 


1. INTRODUCTION AND PRELIMINARIES 


In this paper, we work in the category of smooth objects, so manifolds, maps, 
vector fields, etc, are assumed to be smooth unless it is stated otherwise. 


As is well known, a harmonic map vy : (M,g) — (N,h) between Riemannian 
manifolds is a critical point of the energy functional 


1 
B(y.9)=5 flag? ae. 


The Euler-Lagrange equation is given by the vanishing of the tension filed 7() = 
Trace, Vdy (see [2]), i.e., the map v is harmonic if and only if r(y) = Trace, Vdy = 
0 holds identically. 
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In 1983, J. Eells and L. Lemaire [14] first proposed studying k-polyharmonic 
maps (include biharmonic maps as a special cases). Let us recall the definition of 
biharmonic maps. If a map yọ : (M,g) — (N,h) between Riemannian manifolds 
is a critical point of the bienergy 

1 


Ey. =5f Ito as 


for every compact subset Q of M, then we call it a biharmonic map, where 
T(y) = Trace; Vd is the tension field of y. In 1986, Jiang [19] first computed 
the first variation of the functional to find that y is biharmonic if and only if its 
bitension field vanishes identically, i.e., 


(1) T(y) :— Trace,(V*V* — Vuro) — Trace, R” (dy, T(p))dqg =0, 


where R is the curvature operator of (N,h) defined by 
R'(X,Y)Z = [V¥, VÝ]Z — Vix vj. 


We call a submanifold a biharmonic submanifold if the isometric immersion 
that defines the submanifold is a biharmonic map. By definition, biharmonic 
maps include harmonic maps as special cases, and biharmonic submanifolds can 
be viewed as a generalization of the notion of minimal submanifolds (i.e., har- 
monic (minimal)isometric immersions). Analogously, a Riemannian submersion 
is called a biharmonic Riemannian submersion if the Riemannian submersion is 
a biharmonic map. We call nonharmonic biharmonic maps (respectively, sub- 
manifolds, Riemannian submersion) proper biharmonic maps (respectively, 
submanifolds, Riemannian submersion). 


It is very interesting and important to study the existence problem and classi- 
fication problem which are two fundamental problems in the study of biharmonic 
maps. For the existence problem of biharmonic maps, it would be interesting to 
know if there exists a proper biharmonic map between given two “good” mod- 
el spaces. ( The so-called “good ” spaces include space forms, or more general 
symmetric spaces, or homogeneous spaces, etc.). For the classification problem 
of biharmonic maps, we would especially like to classify all proper biharmonic 
maps between two model spaces where the existence is known. The following two 
classification problems are typical and challenging as follows: 


Chen’s conjecture |12, 13, 11]: any biharmonic submanifold in a Euclidean 
space R” is minimal (i.e., harmonic). 
The generalized Chen's conjecture: any biharmonic submanifold of a Rie- 
mannian manifold (N", h) with non positive curvature Riem" < 0 is harmonic 
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(minimal) (see e.g., [6—13]). 

Though the Chen's conjecture is also known to be true in some other cases (see, 
e.g., [21, 7, 28, 30, 16], it is still open for the general case. For the general- 
ized Chen's conjecture, Ou and Tang (|29]) gave many counter examples in a 
Riemannian manifold of negative curvature. For some recent progress on bihar- 
monic submanifolds, we refer the readers to [1], [4-13], [15-22], [25-31], etc., and the 
references therein. 


Riemannian submersions can be viewed as the dual notion of isometric immer- 
sions (i.e., submanifolds). Based on this, it is interesting to study biharmonicity 
of Riemannian submersions. Biharmonic Riemannian submersions were first s- 
tudied by Oniciuc [22] in 2002. In 2010, Wang and Ou [34] first introduced the 
so-called integrability data and used it to study biharmonicity of a Riemannian 
submersion from a generic 3-manifold, and obtain a complete classification of 
biharmonic Riemannian submersions from a 3-dimensional space form into a sur- 
face. In 2019, Akyol and Ou [1] studied biharmonicity of a general Riemannian 
submersion and obtained biharmonic equations for Riemannian submersions with 
one-dimensional fibers and Riemannian submersions with basic mean curvature 
vector fields of fibers. Moreover, the authors [1] first used the so-called integrabil- 
ity data to study biharmonic Riemannian submersions from (n 4- 1)-dimensional 
spaces with one-dimensional fibers. In 2018, the author [31] studied biharmonicity 
a more general setting of Riemannian submersion with a S! fiber over a compact 
Riemannian manifold. In 2018, the authors in [18] studied generalized harmonic 
morphisms and obtained many examples of biharmonic Riemannian submersions 
which are maps between Riemannian manifolds that pull back local harmonic 
functions to local biharmonic functions. 


Finally, we refer the readers to the following interesting results. In 2011, the 
authors [27] gave complete classifications of constant mean curvature proper bi- 
harmonic surfaces in Thurston's 3-dimensional geometries and BCV 3-spaces, and 
a complete classification of proper biharmonic Hopf cylinders in BCV 3-spaces. 
In 2023, the authors [36] gave a complete classification of biharmonic constant 
mean curvature (CMC) surfaces in 3-dimensional Berger sphere 5?, and a com- 
plete classification of proper biharmonic Hopf cylinders in S2. On the other hand, 
in 2010, the authors [34] gave complete classifications of Riemannian submersions 
from a 3-dimensional space form. In 2023, the authors [35] classified all proper 
biharmonic Riemannian submersions from BCV 3-diemnsional spaces into a sur- 
face, and proved that a proper biharmonic Riemannian submersions from a BCV 
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3-space exists only in H?x R — R?, or, SL(2, R) — R°. In a recent paper [36], the 
authors also gave complete classifications of biharmonic Riemannian submersions 
from a 3-dimensional Berger sphere and proved that a biharrmonic Riemannian 
submersion from a 3 Berger 3-sphere 5? into a surface has to be harmonic. 


In this paper, we study biharmonic isometric immersions of a surface into 


M? x R and biharmonic Riemannian submersions from product spaces M? x R 


into a surface. We give a complete classification of proper biharmonic isometric 


immersions of a surface with constant mean curvature into product spaces M? x R 
. More precisely, we prove that a proper biharmonic surface with constant mean 


curvature H in M? x R is a part of S! (zziz) x R in S? (an) x R. We also 


obtain a complete classification of proper biharmonic Hopf cylinders in M? x R. 


On the other hand, we give a complete classification of biharmonic (including 


harmonic) Riemannians submersions 7 : M? x R — (N?,h) from M? x R toa 
surface and construct many family of proper biharmonic Riemannian submersions 
T: M? xR >R. 


2. BIHARMONIC ISOMETRIC IMMERSIONS OF A SURFACE WITH CONSTANT 
MEAN CURVATURE SURFACES INTO M? x R 


Biharmonic surfaces in 3-dimensional space forms have been completely clas- 
sified ([20], [11], [9], [10]). In [27], the authors classified all proper biharmonic 
CMC surfaces in 3-dimensional BCV spaces and Thurston’s 3-dimensional ge- 


ometries. Biharmonic CMC surfaces in M?(c) x R, where M?(c) denotes a space 
form, has also been studied in [15] and [17]. In the recent paper [36], the authors 
give a complete classification of biharmonic constant mean curvature surfaces 
in 3-dimensional Berger sphere S?. In this section, we will study biharmonic 
isometric immersions of a surface with constant mean curvature surfaces into 


M? x R and proper biharmonic Hopf cylinders in M? x R, where M? is a general 
2-dimensional manifold. 


It is well known that a 2-dimensional Riemannian manifold (M?, g) can be 
locally expressed as (R°, e??)da? + dy?) with respect to local coordinates (z, y) 


in R?. Without loss of generality, locally, we identify product spaces M? x R with 
(IR?, e?» (o dz? + dy? + dz?) with respect to local coordinates (x,y,z) in R. 
We adopt the following notation and sign convention for Riemannian curvature 


operator: 


(2) RQOGY)Z e Vx V9Z — Vy Vy4 — VisnZ, 
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and the Riemannian and the Ricci curvatures: 
R(X, Y, Z,W) = g(R(Z,W)Y, X), 
3 


(3) Ric(X, Y) = Trace Ho, BOY ei, X, ei) = UIX e;)e;, Y). 
i=l 


3 
=] 


l 


It is not difficult to check that (E, = e Play) 9. E = i Ez = at form an 
orthonomal frame on M? x R = (R? x R, e??) dx? + dy? + dz?). 


With respect to this frame , a straightforward computation shows that 


(4) |E, Eo] = pÉ, all other |E; E;] = 0, i, J = 1,2,3. 


Let V denote the Levi-Civita connection of M? x R, then we can check that 


(5) Vg, Ei = —py E», V pg, E = ae all other Vg, Ej = 0, i, J = 1, 2,3. 
A further computation gives the possible nonzero components of the curvatures: 
(6) Ria = g( (Ei, Ey) E2, E1) = —Pyy — B 

all other Fou = g 2a 2) Ey Ei) = 0, i,j, k,l = 1,2,3. 


and the Ricci curvature: 


(7) 
Ric (£4, £1) = Ric( E2, E2) = —Pyy — på, all other Ric (Ej, Ej) = 0,2, j = 1,2,3. 


The following equation for biharmonic hypersurfaces in a generic Riemannian 
manifold will be used in our study of biharmonic surface in the product space. 


Theorem 2.1. ({26]) Let o : M" —> N""*! be an isometric immersion of 
codimension-one with mean curvature vector n = HE. Then q is biharmonic if 
and only if: 
(8) AH — H|AP + HRic" (£, £) = 0, 
2A (grad H) + "grad H? — 2 H (Ric" (£))' = 0, 

where Ric : T,N —> T,N denotes the Ricci operator of the ambient space 
defined by (Ric (Z),W) = Ric "(Z, W) and A is the shape operator of the hy- 
persurface with respect to the unit normal vector £. 


We now study biharmonic surfaces with constant mean curvature H in the 


product spaces M? x R. 


Theorem 2.2. A surface with constant mean curvature H in the product space 


M? x R is proper biharmonic if and only if it is a part of S! (zzz) x R in 


g? (an) x R. 
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3 3 
Proof. Let (e = 2 at Ei, €2 = » ak, E= 2 al E;) be an adapted orthonor- 


mal frame with £ T normal Xp the ses. By a simple computation, we 
use the Ricci curvature (7) to obtain Ric (€,€) = R(1—(a3)?), (Ric (€))' = 
—Ra3(ajei + a3e2), where R = —py, — p;. From these and using biharmonic 
surface equation (8), one finds that a surface with constant mean curvature H is 
biharmonic if and only if 


(9) —H ||A|? — R(1— (a3))] = 0, Raja$H = 0, and Ra3a3H = 0. 
This has solution H = 0 implying that the surface is minimal, or, 
(10) JAP = R(1 — (a), Rata = 0, and Rates = 0, 


where R = —py-— in. 

We solve Eq. (10) by considering the following cases: 

Case I: R = —py, — p; = 0. In this case, we have |A|? = 0 and hence H = 0 
meaning that the surface is minimal. 

Case II: R = —pyy — p, # 0. In this case, by the last two equations of (10), we 
have either a3 = 0 or a? = a3 = 0. This is discussed by the following two cases. 


For Case I-A: a3 = 0, applying the first equation of (10) we conclude that 
(11) |A|? = 


Clearly, a3 = 0 means that the normal vector field of the surface X is of the 
form € = al E; + a2Es. This implies that the normal vector field € is always 
orthogonal to E3 = 2. Consequently, we can choose an another orthonormal 
frame {e; = aEı + bEo,e2 = E3,€ = bE, — aE} adapted to the surface with 
a? + b? = 1 and € being the unit normal vector filed. By a straightforward 


computation using (5) we have 


(12) Ve = (aei(b) — bes (a) — apy e1, Ve, = (aea(b) — bea(a)) e1. 
With respect to the chosen adapted orthonormal frame, a further calculation 
gives the second fundamental form of the surface X as 

(13) 

h(ei, e) = —(Va£, €1) = —{ae,(b) d bei(a ) m apy}, h(ei, €2 E —(V i; €2) = 0, 
h(€2, é1) = xi €1) mg (ae2(b) m bes(a )) ’ h(e2, e2) = =A, ae e2) = = 0. 

It follows from (13), the symmetry h(e1,e3) = h(es,e1) = 0 that e2(b) = ez(a) = 

0, h(e1,61) = 2H, and hence |A|? = 4H?. From this and using (11), we have 


(14) R(E, E», E E2) = —Pyy — D, eq eH. 


This implies that M? has constant Gauss curvature KM = 4H? > 0. It fol- 
lows from a well-known fact in the differential geometry of surfaces that M? 
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can be identified with 2-sphere S? (zin 


c = AH?. Therefore, it follows from Theorem 2.1 in [27] that a constant mean 


) with positive constant Gauss curvature 


curvature surface in S? (an) x IR is proper biharmonic if and only if it is a part 


of S! ( x R in S? (stn) x R. 


zm) 


For Case II-B: a? = a3 = 0 and hence a$ = +1. It follows that Span{e1, e2} = 
Span{ £,, Ej). This means that the surface is M? in M? x R, which is a minimal 
isometric immersion . 


From which, the theorem follows. 


Theorem 2.3. Let * : M? x R > M? be a Riemannian submersion with 
n(x,y,z) = (x, y) defined by the projection onto its first factor and B : I > M? 
be an immersed regular curve parametrized by arc length. Then the Hopf cylinder 
X = Uer '(8(s)) is a proper biharmonic surface in M? x R if and only if it is 
a part of S! ( 


ME E 2 . | 
m) x RR ins zi) x IR with nonzero constant mean curva 


ture H. Furthermore, this holds implying that M? x R has to be S? (stn) x 


A 


? 


i.e., M? identify with a 2-sphere S? (zin) with positive constant Gauss curvature 
ped 


Proof. Consider the Riemannian submersion m : M? x R = (R?, ere qg? + 
dy? + dz?) — M? = (R?, cP) da? + dy?) with m(z,y, z) = (zx, y) and let 8: I > 
M? = (R?, e??@) dx? + dy), B(s) = (x(s), y(s)) be an immersed regular curve 
parametrized by arc length with the geodesic curvature k,. We take the horizontal 
lifts of the tangent and the principal normal vectors of the curve 6 : X = aE4--bEs 
and £ = bE, — aE, (where a = e"9)5/, b = y' and a? + b? = 1) together with 
V = E; to be an orthonormal frame adapted to the Hopf cylinder (see also [26]). 
A direct computation using (7) gives: 


(15) Ric (E, £) =k= ~Pyy 7 97; Ric (£, X) = 0, Ric (£, V) = 
We can check that the the geodesic torsion of the lifting curve 7~'(3(s)) 


(16) Tg = —(VxV,£) = —(Vag, dog, E3, bE; — a ES) = 0. 
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It follows from Eq. (16) in [26] that the surface X = Uierr !(8(s)) in M? x Ris 
biharmonic if and only if 


kj — kg(k? + 272) + kj Ric(£, €) = 0, 

3k! ky — 2k, Ric(£, X) = 0, 

k Tg + kRic(6, V) = 0. 
We substitute (15) and (16) into the above equation to get 
(17) kj —ki+k,R — 0, and 3k; k, = 0, 
where R = —py, — B. 
Solving (17) we have k, = 0 implying that the surface X is minimal surface, or 8 
has constant geodesic curvature k, and k? = R = —Pyy— Py > 0. We use the fact 
from [26] (Page 229) to conclude that the mean curvature of the Hopf cylinder 
given by H — "i and hence |A|? = k? + 277 = k = R = -py = p2. From these 
we conclude that the Hopf cylinder X = Usern !(0(s)) is proper biharmonic if 


and only if 
(18) 
= pd: 
H? = È = ^. = constant > 0, |A|? = R= -pyy — p, = constant > 0. 


We apply the characterizations of Hopf cylinders in M? x R given in Theorem 
2.2 to obtain the Theorem. 


Corollary 2.4. A totally umbilical surface in M? x R is biharmonic if and only 


if it is minimal. 

Proof. Since a totally umbilical biharmonic surface in 3-dimensional Riemannian 
manifolds must have constant mean curvature H ({27]), we have from Theorem 
2.2 to conclude the only potential totally umbilical proper biharmonic surface is 


a part of S! ( x R in S? (zin) x R. Clearly, it is not totally umbilical. 


eod - 
2/2|H| 


3. HARMONIC AND BIHARMONIC RIEMANNIAN SUBMERSIONS FROM M? x R 


In this section, we study biharmonicity (including harmonicity) of submersions 
from M? x R — (N?, h) and also derive a classification of biharmonic (including 
harmonic) Riemannian submersions from M? x R by identifying M? x R with 
(IR?, e?» (o: dx? + dy? + dz?) in local coordinates. 


Let 7 : M? x R > (N?,h) be a Riemannian submersion from M? x R to 
a surface with an orthonormal frame {e1, €», e3} and e3 being vertical. By a 
computation similar to Remark 1 in [35], we have the following (19)-(25) (see 
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[35] for details). 
The Lie brackets given by 
(19) [e, e3] = faea + Kies, [e2, 63] = — fae + K2e3, [e1, 02] = fier + foes — 2063. 


where (fi, fo, f3,41, K2, c] is the (generalized ) integrability data. Moreover, 
f = 0 if and only if the above frame is adapted to 7. 
The Levi-Civita connection can be given by 


(20) 
Vaer = —fiez, Veal = fier — 063, Ve,e3 = 065, 
Ve.€1 = — f2€2 + 023, Ves€2= foe, Vez€3 = —7€1, 
Ve,€1 = —Ki€3 + (0 — fz)e2, Vege2 = —(o — fz)€1 — K2e3, Vezg€3 = K1€1 + Ka€2, 


3 
We denote by e; = >> a2E;, i = 1,2,3 and use (5), (6) and (20) to check that 
j=l 
the Jacobi identities as 


(21) e3( fi) + (Ki + fo) fs — e1( fs) = 0, e3(f2) + (Ka — fi) fs — e2( f3) = 0 
2e3(0) + Ki fi + Rafa + es(K1) — ex(&2) = 0, 


and the terms of the curvature tension as follows 


RV (e, €3, €1, €2) = —ei(0) + 2k10 = —a3a3R, 
RM (ei, €3, €1, €3) = €1(K1) +0? — K? + kofi = (a3)? R, 
R™ (e, €3, €2, €3) = €1(K2) — e3(0) — Ki fi — Kio = —aja3R, 
(22) 4 R” (ei, e2, €1, €2) = ex(f2) + ex (fi) — fÈ — f2 + 2f30 — 3o? = (a3) R, 
RM (e, 65, 65, 63) = —€2(0) + 2K20 = aja$R, 
RM (eg, €3, €1, €3) = €2(K1) + ea(o) + Ko f2 — Kiko = —aja3R, 
RV (e2, es, €2, €3) = a? + €2(K2) — Ki f2 — KZ = (at)? R, 


where R = —py, — B 
We have Gauss curvature of the base space as follows 


(23) KY = e(f2) — ea(fi) — f? — f2  2fao. 
Moreover, 
(24) e3(K™) = es{e1(fo) — ea( f1) — f? — f2 + 2fso} — 0. 


When f3 = 0, Gauss curvature of the base space becomes 


(25) K™ = e(fg) — efi) - fi — fa. 


10 ZE-PING WANG* AND YE-LIN OU** 


3.1. Harmonic Riemannian submersions from M? x R. It is a well known 
fact that harmonic maps are always biharmonic maps. Therefore, it would be 


interesting to study harmonicity of Riemannian submersions from M? x R be- 
fore studying biharmonicity of them. In [37], the authors classified all harmonic 
Riemannian submersions from Thuston's 3-spaces, BCV 3-spaces and Berger 3- 
sphere. In this subsection, we will study harmonic Riemannian submersions from 
M? x R. In next subsection, we we will classify all proper biharmonic Riemannian 


submersions from M? x R. 


Now we are going to give the following classification of harmonic Riemannian 


submersions from M? x R. 


Theorem 3.1. A Riemannian submersion x : M? x R > (N?,h) is harmonic if 
and only if it is the projection 7 : M? x R > M?. 


Proof. Let m : M? x R — (N?,h) be a Riemannian submersion with an or- 
thonormal frame (ei, e2, ea), e3 being vertical, and the (generalized) integra- 
3 

bility data (fi, fo, fs, K1, 2,0]. Denoting by e; = MjajE;, i = 1,2,3. It fol- 
j=l 

lows from Proposition 2.2 in [37] and (22) that if a Riemannian submersion 
7: M? x R = (R^ x R, 7?) dx? + dy? + dz”) — (N?,h) is harmonic, then 
we have 


(26) KN = 30? + (a3)?R, 


where R= —py, = p? and KY = e;( fz) + ex(i) — f? — f2 + 2fso. 


Since Gauss curvature of M? equals R = —py, — i it follows that if R = 


—Pyy — p, = 0 then M? is a Euclidean space R7, i.e., the product space M? x R 
is actually a Euclidean space R? x R = R?. On the other hand, we use the 2nd 
equation and the 4th equation of (26) to conclude that o = 0 and hence Gauss 


curvature of the base space K^ = 0, which implies that the base space N? can 


be viewed as a 2-dimensional Euclidean space R?. In the case, a harmonic Rie- 


mannian submersion 7 : M? x R — (N?,h) from a product space M? x R exists 
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only in the orthogonal projections 7 : IR? = R? x R > R°. 


From now on, suppose R = —py, — Py z 0. In this case, we first use the 2nd 
equation and the 7th equation of (26) to have (aĵ)? = (a3)?. On the other hand, 
using the 3rd equation and the 6th equation of (26) we have a?a3 = 0. From these, 
we have a? = a3 = 0 and hence (a3)? = 1. We apply the 2nd and the 4th equation 
of (26) to get c? = 0 and hence Gauss curvature of the base space KY = R equals 
Gauss curvature of M?. This implies that N? is locally isometric to M?. Since 
(a3)? = 1, one finds that the orthonormal frame (e; = Ej, e3 = E», e3 = Es], 
together with (5), is an adapted frame to the Riemannian submersion m with 
e3 = E; being vertical. Combining these, a harmonic Riemannian submersion 
na: M? x R > (N?,h) is actually the projection 7: M? x R  M?. 
From which the theorem follows. 


3.2. Biharmonic Riemannian submersions from the product spaces M?x 
IR. In this subsection, we will classify all proper biharmonic Riemannian submer- 


sions from M? x R, which are not harmonic. 
The following lemma was obtained in [34] which will be later used in the rest of 
the paper. 


Lemma 3.2. ([34]) Let  : (M?,g) —> (N?,h) be a Riemannian submersion 
with the adapted frame (ei, es, e3} and the integrability data fi, fa, K1, K2 and a. 
Then, the Riemannian submersion m is biharmonic if and only if 

(27) 


= AM gs — 25 fie o) — Ka » (ei (fi) — fifi) + ra (-*" as n) ub 
-AAMk 4.2 > fiei(ki) + à Y (ef) — Ki fa) + Ko (-i* + 2 7) ed). 


i=l i=l i=l 
where KN = Rẹ asor = ei(fa)—-ex( fa) - f2—f2 is Gauss curvature of Riemannian 
manifold (N?, h). 


The following proposition was derived in [35] 


Proposition 3.3. (see [35]) Let x : (M?,g) — (N?, h) be a Riemannian submer- 
sion from 3-manifolds with an orthonormal frame (e1, €2, e3) and ea being vertical. 
If Vee = 0, then either V.,€2 = 0; or Ve,e2 Æ 0, and the frame {e1, €2,e3} is 
adapted to the Riemannian submersion m. 


We apply Proposition 3.3 to have the following theorem which will be used to 
prove our main result 
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Theorem 3.4. Let 7: M? x R = (RË, e?» 9)dz? + dy? + dz?) — (N?,h) be a 
Riemannian submersion from product spaces with R = —pyy — p, z 0. Then, 
there exists such an adapted frame (e, = al E, +a? En, e2, e3} of the Riemannian 
submersion t with ea being vertical. Moreover, if E is not vertical, then V.,€2 4 


0, i.e., fo X 0. 


Proof. Clearly, if Es is tangent to the fiber of the Riemannian submersion 7, i.e., 
e3 is perpendicular to both E, E», then any basic vector field can be expressed 
as the form e = a?E, + b? E», and a? + b? = 1. 


From now on, we need only to suppose that E; is not vertical, i.e., e3 is not 
parallel to E3. Then, the vector filed ey = e3 x Es; is horizontal and hence 
(e1, £3) = 0. From this, one can obtain a defined orthonormal frame (ei, e = 


3 
e3 X ej, eas) on M?. If denoting by e; = >> alE;i = 1,2,3, together with 
j=l 


(e1, £3) = 0, then the vector horizontal filed e} can be expressed as the form 
e1 = al E, + a? E and hence a? = 0, (at)? + (a2)? = 1. From these, we have the 
following 


(28) a = 0, a3» +1 and a3 Æ 0. 
Moreover, one has also the following equalities as 
(29) fi = 0, Wa = 0. 


In fact, by a straightforward computation we get 


3 3 MN 
(30) Veer = Ved) Ei) = } e(a1) Ei 9, aja Vg, Ei. 
i=l i=1 ij-l 
In addition, using(20), the above has another expression as 
3 
(31) Veer = —fie2 = -fı S dE. 
i=1 
Equating (30) and (31) and making a comparison of the coefficient of Es, we 
obatin 
3 . 
-= fia} EE (—-R b» aj Ei, E3) = (Vei; E3) 
(32) 3 e 
= (X e(a)E * Yo aiai V g, Ei, E3) = ex(a1) = 0, 
i-l ij-l 


which has been used (5) and aj = 0. This deduces fi = 0 for a3 Æ 0, from which 
we get (29). 
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Applying (5), (20) and a? = fı = 0 and performing a computation in the same 
way to those used computing (30)-(32), we have 


e1(a3) = —eaj, 
e1(a3) = 043, 
e>(a3) = 0, 

(33) ex(a3) = 0, 
ea(a3) = —K2a3, 
e3(a3) = K203, 


Since Ve e1 = 0, we apply Lemma 3.3 to have either Ve,€2 # 0, and the frame 
1€1,€2; 63) is adapted to the Riemannian submersion 7, or V.,€2 = 0. For the 
latter case: V.,€2 = 0, i.e., f? = 0, we will show that if a3 = 0, then the chosen 
frame {€1, €2, e3} is also adapted to the Riemannian submersion m whilst the oth- 
er case a3 # 0, +1 and f? = 0 can not happen. 


In fact, first of all, we show that the case a3 # 0,+1 and f; = 0 can not 
happen. If otherwise, substituting f; = 0 into the 8th equation of (33), together 
with a3 Æ 0 and a3 Æ 0, +1, we get o = 0. Furthermore, substituting these and 
fi = fo = 0 into the 4th equation of (22), we have R = —Pyy— Py = 0 contradicting 
the assumption R # 0. Secondly, we now consider the case fz = aj = 0. Since 
a = aj = 0, it is very easy to see that (a3)? equals 1. We apply the 2nd equation 
and the 7th equation of (33) separately to obtain o = 0 and hence f3 = 0. 
Combining these, for the case a3 = 0 and aj = fı = f» = 0, the frame {e1, e2, e3} 
is an adapted frame of the Riemannian submersion. 
From which, we obtain the theorem. 


Remark 1. Let m : M? x R = (R?, e? (da? + dy? + dz?) — (N?, h) be a Rie- 
mannian submersion. Then, we have the following cases: 


(a) E R mec p? = 0, then the potential product space is IR^. We ap- 
ply Theorem 3.3 in [34] to conclude that biharmonic Riemannian submersion 
T : R? — (N?, h) has to be harmonic. 


(b) If (a3)? = 1, i.e., E3 is vertical, then we use (5) to obtain the tension of 
the Riemannian submersion 7(7) vanishes, i.e., 7 is harmonic. 
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Remark 2. Let n : M? x R = (R?, e??@¥)dx? + dy? + dz?) — (N?, h) be a Rie- 
mannian submersion with e3 being vertical. Suppose that Æ; is not vertical (i.e., 
a; # +1) and R = —py, — p; # 0. From the proof of Theorem 3.4, we can 
choose such an orthonormal frame (e; = aj E; + a? En, €2,€3} on M? x R with 
aj = f, — 0, a$ Z +1 and a3 Z 0. Applying Theorem 3.4 and together with the 
proof of Theorem 3.4, one can further find the following: 


(i) fa #0. In this case, the chosen frame fei = al Ei + a Es, e2, e3} is actually 
adapted to the Riemannian submersion 7 with a? = fi = f3 = 0. 


(ii) the case fo = 0 and a3 £ 0, +1 can not happen. 


(iii) fo = 0 and a} = 0 and hence aj = +1. We show that Riemannian sub- 


mersions 7 must be M? x R — R?. In fact, note first that the chosen frame 
, 


{e, = aLE, + af? E>, e2,€3} is adapted to the Riemannian submersion m with 
a}? = fi = f = 0. For aè = @} = fi = fə = fs = 0 and a} = +1, by a 
simple computation using the 3rd, the 6th and the 7th equation of (33) gives 


o = ko = 0. Moreover, from these, we use the 2th equation of (22) to conclude 
kı # 0. In addition, it is easy to find from (23) that Gauss curvature of the base 
space KN = ei(fz) — ex(fi) — f? — f2 + 2fs0 = 0 and hence N? is a Euclidean 
space R*. Therefore, Riemannian submersions 7 must be M? x R > R?. 


(iv) On the condition of (iii), there exist many family of proper biharmonic 
Riemannian submersions 7 : M? x R — R°. Moreover, with respect to local 


y ; : ; i : 2 
coordinates, a proper biharmonic Riemannian submersion 7 : M? x R > R° can 
be locally expressed as 


a: M? x R = (R3, e» 9q52 + dy? + dz”) > (R?, dy? + dz”), m(v,y,z) = (y,z), 


and the function p(z, y) satisfies PDE (90) and p, ¥ 0. 


In fact, it is observed that the orthonomal frame (ei = alE, + af E>, ej = 


Es, e3 = alE, + alE3] on M? x R is an adapted frame to the Riemannian 


submersion 7 with e3 being vertical. Note also that the basic vector field e; = 
al E +a% E; and the vertical field e3 = al E; 4-a2Es are tangent to M? and another 
basic the vector field e3 = Es is perpendicular to M? (i.e., e = Es parallel to R). 


Since V,,e, = 0, it is easy to see that an arbitrary integral curve of e, is geodesics 


on M? c M? x R. It follows from a well-known fact in the differential geometry 
of surfaces that one can parametrize M? by r — r(u,v) so that the v—curves 
are the integral curves of e; and the u—curves are the orthogonal trajectories 
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of the v—curves, and the Riemannian metric on M? can be also expressed as 
g = e24)du? + dv? by reparametrization with e3 = e-*" 2, e; = 2 and 


hence R = —quy — q = —Pyy — p, on M?. It turns out that with respect to local 


coordinates (u, v, z), the Riemannian submersion m : M? x R > R° can be locally 
expressed as 


T : M? x R = (R3, 0») du? + dv? + dz”) > (R?, dv? + dz”), n(u, v, z) = (v, z). 


Since the Coordinate transformation (x,y) — (u,v) is an isometric transforma- 
tion M? — M? and it is a convention to use coordinates (x,y) in place of (u,v), 
the above Riemannian submersion 7 : M? x R > R? can be also locally expressed 
as 


T M? xR= ( R3, ero) dg? + dy? RE dz?) xc ( R2, dy? + dz’), n(z, y, z) m (y, 2). 


In addition, by Proposition 3.7, Remark 3 and Example 1, one sees that there ex- 


ist many family of proper biharmonic Riemannian submersions 7 : M? xR > R’, 
where M? x R satisfies certain conditions determined by p(x, y) solving (90). 


We now give a characterization of Riemannian submersion m : M? x R — 
2 
(N*,h). 


Theorem 3.5. If 7: M? x R > (N?,h) is a proper biharmonic Riemannian 
submersion from the product space, then 
(a) The target surface is flat, and locally the map can be expressed as 


T: M? x R 2 (R3, e?» 9q52 + dy? + dz”) > (R?, dy? + dz”), r(x, y, z) = (y,2), 
where p(x, y) is a function solving PDE (90) and p, 4 0, or, 

(b) With respect to another local coordinates, it is a Riemannian submersion 
1: M? x R= (IR?, 22M dy? + dv? + dz?) > (R?, dv? + ede’), 

nlu, v, 2) = (wv, 0) = (v, FAW + Az + B), 


where Q(v) = In|tana(v)|, A = In|sina(v)| — In A depend on only v , B and 
A> 0 are constants, and a nonconstant function a(v) solves the following ODE 


(34) 


(35) a” sin a cos? a + cos o(sin? a + 3)oo" + sin a(2 cos? a + 3)o/? = 0. 


Proof. For the proof of the theorem, we will use a local expression of M? x R with 
the form (IR? x R, e??) dx? +dy?+dz?) in local coordinates, where M? has such a 
local expression the form (R°, e??(*)dx?+dy?). Let V denote the Levi-Civita con- 


3 
nection on M? x R = (IR?, c??@) dx? + dy? + dz’) and by e; = > a! Ej, i = 1,2,3. 
jal 


For considering proper biharmonic Riemannian submersions y : M? x R > 
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(N?,h), by Remark 1, we need only to consider R = —py, — Py = 00 
and aj # +1 in the rest of the proof of the theorem. Therefore, by Theo- 
rem 3.4, there is an adapted frame (ei, e», e3} to the Riemannian submersion 
m with e; = alE, + a?Ez, e3 being vertical, a? = fı = f3 = 0, and hence 


e3( fi) = e3(f2) = 0 (see [34]). 


By Remark 2, we now just need to consider the following two cases: 
Case I: ał = 0. In this case, we have a3 = +1 since af = 0. Noting that f3 = 0, 
we use the 8th equations of (33) to have fy = 0. By (iii) and (iv) of Remark 2, 
we immediately give the statement (a). 


Case II: a3 Z 0, X1. This means a3 4 0, +1 since a? = 0. In this case, we show 
that for f; # 0, there exist also proper biharmonic Riemannian submersions from 
certain product spaces M? xR. 

By Remark 2 and combining the above, we have the following conditions: 


(36) as, 05 7 0,1, f2 #0, aj = fi = fs =, ex (JU) = e3( f2) = 0. 


This case corresponds to the statement (b) which is obtained by the following two 
steps: 


Step 1: show that e3(f2) = e»s(&1) = es(&1) = €2(0) = e3(0) = e(R) = 
e3(R) = ka = 0, «1o £0. 


Firstly, we have o # 0 by the 8th equation of (33). 


Secondly, show that 4&9 = 0 and hence e2( f2) = ea(f2) = ex(e) = e3(0) = 
€9(K1) = €3(K1) = e3(R) = e2(R) = 0). 


We use the 5th and the 6th equation of (33), (24) and (36) and apply es to 
both sides of the 8th equation of (33) to get 
(37) — fok2a3 = ea(c)a$ + 63a. 


One substitutes the 8th equation of (33) into (37), together with (36), and sim- 
plifies the resulting equation to have 


(38) €3(0) = —Kaf2 — KK. 
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Substituting the above result into the 6rd equation and the 3th equation of (22) 
separately, together with (36), we have 


(39) €2(K1) = 2K1K2, e1(&2) = — Kə f2. 


A straightforward computation using the 7th, the 8th equation of (33), the 5th 
equation of (22), (37) and (36) gives 


1 
(40) e3(0)a3a3 = — k0 = =z), (fa + K1)a3a5 =o. 


Together with the 8th equation of (33) and (36), we multiply the left and the right 
hand sides of the 7th equation of (33) by f2a3 and ca$ respectively, to obtain 


(41) a = Kı fa. 
Substituting this into the 7th equation of (22) with a} = 0 yields 
(42) €9(K) = Kå. 


We use the 6th equation of (33) and (24), and apply e3 to both sides of the 4th 
equation of (22) to get 


(43) —Ó6cea(c) = 2K,a3a3R + (a3)? ea( R). 
On the other hand, a direct computations using (20), (22), (33) and (36) gives 


V R(e1,€2, 61, 62; 63) = Ve, R(e1, €2, €1, €2) 
(44) = eg(R(e1, €2)€2, e1) — (R(e1, e2)e2, Ve,e1) — (R(e1, e2)V «462, €1) 
—(R(ei, Ves€2)€2, €1) — (R(V 01, €2)€2, €1) = (a$)?ea(R). 


In the same way, we can check the following 

(45) V R(e1, €2, €3, €1; €2) = ajaSea(R), V R(e, €2, e2, 63; €1) = 0. 

Using (44) and (45) and applying the Second Bianchi-identity, we have 

0 = V R(e1, €2, €1, €2; €3) + V R(e1, €2, 63, 61; €2) + V R(e1, €2, €2, €3; €1) 


46 

M6) — (a2)e(R) + adade(R), 

that is, 

(47) (a3)e3(R) + adades(R) = 0, 


where R = —pyy — py # 0. 


We use the 8th equation of (33), (36) and (40), and apply e3 to both sides of 
the 1st equation of (41) and simplify the resulting equation to have 


(48) es(&1)(a3)? = —2k0, €3(K1)a3 = 2es(c)a3. 
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Applying es to both sides of the 2nd equation of (22), a straightforward calculation 
using the 3th equation of (33), (36), the 1st equation of (39) and the 5th equation 
of (22) gives 
(49) 

€2€1(K1) = 2€9(K1) — 20e2(0) + (a3)? e3(R) = 4213 — 2&0? + (a3)? es( R). 


We apply ei to both sides of the 1st equation of (39) and use the 2nd equation 
of (22), (36) and the 2nd equation of (39) to give 


(50) €1€2(K1) = 2K1€1(K2) + 2K 2€1(K1) = —4ko0? + 2K? K2 + 2ka(a3 ^ R. 


Noting that foe? — 20e3 = |e1, e2] = e1€2 — €se1, we can compute using (49) 
and (50) to obtain 


pu (51) 
= foe2(K1) — 20€3(K1) = €1€2(K1) — e2€1(K1) = —2&20? — 2K? Ky + 2K2(a3)?R — (a3)? e3( R). 


Substitute (39) and (41) into (51) to have 


(52) (a3)*e(R) = —4k go? — 2K? Ky + 2k2(a3)?R + 20ea(14). 

— A further computation using the Ist equation of (48) yields 

(53) 

~ (a3)(a3)?e2(R) = —4k20°(a3)” — 2r1K2(03)” + 2&2(a5) (a3) R + 20e3(%1)(a3)” 
pir = —Ak36? (a3)? — 2&1&»(a3)? + 2&3(a3)? (a)? R — A&30?. 


Substituting (43) and (53) into (47) to compute, we obtain 


0 = aga (a3a$e? (Ft) + (a3)^es(R)) = (a3) (a3)*e3() + (a3)? ajases(H) 


= —4&30? (a3)? — 2&2&»(a3)? — 6oea(c)a3a3. 


= (54) 
Substitute the 1st equation of (40) and the 7th equation of (33) into the above 
equation together with (a3)? + (a3)? = 1 and fs = 0 to obtain 


(55) 290" (a3)^ = 0, 

this implies kə = 0 since a3 # 0 and ø zz 0. We substitute x = 0 into (38), 
(39),(40), (43), (47) and (48) separately to have ky = e2(7) = e3(0) = €2(K1) = 
€3(K1) = e3(R) = ex(R) = 0. Applying ez to both sides of the 8th equation of 
(33), we immediately obtain e3( f2) = 0. 


Thirdly, we show that M? x R can be also represented as (IR?, e?9 qj? - dv? + 
dz?) with respect to another local coordinates and kı = K1(v) = — sin? a(v)Q'(v), f; = 
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fo(v) = — cos? aQ', o = a(v) = — sina cos aQ', o/(v) = —e = sina cosaQ'(v), and Q(v) = 
In | tan a(v)]. 
Since aj — 0, we can assume that 


€1 = cos 0E, + sin OE», 
(56) e> = —sin@ cosa E4 + cos 0 cos a Es + sin a E3, 


e3 = sin f sin aE — cos 0 sin a Es + cos a E3, 


where the functions a,@ : M? x R > R, denote the angle between e3 and Es, 

between e4 and £i, respectively. 

Denoting by E, = sin QE; — cos0E», E; = cos0 E, + sin& E», then (56) can be 
T rewritten as 


Y (57) e = Ey, e = —cosak, +sinaks3, e3 = sina E, + cosaks. 


Noting that sina = a3,cosa = a3, we use (33) and k2 = 0 to have e2(sin a) = 
e3(sina) = 0, e2(cosa) = e3(cosa) = 0 implying e2(a) = e3(a) = 0. A straight- 
forward computation using (5), (20) and f3 = 0 gives 


K1€1 = Ves = Ve, (sin 0 sin aL; — cos @ sin qE + cos Ei) 
mas = e3(sin 0 sina) E; — ea(cos 0 sin a) E» + e3(cos a) Es 
iE + sin 0 sin a Ve, E, — cos0 sin aV e, E» + cos aV, E3 
= = cos Üsin we3(0) E, + sin 8 sin aea(0) Es 
(58)  -F(sin 8 sin a)?V p, E, — sin 0 cos 0 sin? aV p, E, + sind sina cos aV p, Ey 
— sin 0 cos 0 sin? aV p, E + (cos 0 sin a)? V g, E; — cos0 sin a cos aV p, E» 
— + sin Ó cos asin oV p, £3 — cosÓ cos a sin aV p, E3 + cos? aV g, Es 
e = sina(es(0) — sin 0 sin ap,){cos GE, + sin 0 Es 
= sina(ea(0) — sin 0 sin o,)es, 


which implies 
(59) kı = sina (e3(0) — sin 0 sin apy) . 


Similarly, we can check the following 
(60) 
—kK1€3 + oez = Ve; = Veg (sin 0E + cos 0E2) = (e3(0) — sin 0 sin apy ){— sin OE; + cos 0E»), 


and 
(61) 
— foeg + aes = Vei = Ve (sin OE; + cos 0 E2) = (e2(0) + sin cos ap,){— sin 0E, + cos 0E}. 
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Substituting (56) into (60) and (61) separately, we get 


(62) o = cosa (es(0) — sin 0 sin apy) , 
(63) o = — sina (e2(0) + sin 0 cos apy), 
and 

(64) fo = — cos a (e»(0) + sin 8 cos apy) . 
We compare Eq.(62) with Eq.(63) to conclude the following 
(65) cos aea(0) + sin ae2(0) = 0, 
which implies 

(66) Ex(0) = Ẹ(0) — 0, 


it follows that the function 0 doesn't depend on the variable z and hence (£4 = 
sin GE, — cos 0 E2, E; = cos 0E, + sin 0E = ej) is an orthonormal frame on M?. 
Since Vg, E; = Vee = 0, then an arbitrary integral curve of e4 is geodesics 
on the surface M? C M? x R. It follows from a well-known fact in the dif- 
ferential geometry of surfaces that we can reparametrize M? by r — r(u,v) so 


that the v—curves are the integral curves of e1 and the u—curves are the or- 
thogonal trajectories of the v—curves, and the Riemannian metric on M? can 


be locally represented as g = e? du? + dv? with E, = ec 2. By = e, = 
2 with respect to local coordinates (u,v). Combining these, a product space 


M? x R can be also locally expressed as (RÌ, e?9 9 qu? + dv? + dz?), where 
u = u(r,y),v = v(z,y). It is observed from (57) that the orthonormal frame 


fe, = Eo, e3 = — cosa E, + sinaEs, e; = sin a F; + cosa Es} is adapted to c with 
e3 being vertical. A direct computation gives 
Va E; = —qy Ey, V g, E» — qu EA, 


(67) all other Vg,E; 20, i,j = 1,2,3, 


where E4 denotes E3. 
Clearly, 
(68) 
—qw — & = 9(R(E1, E) Ey, E1) = R = g(R(Ei, Ex) Ey, E1) = —pyy — v5. 


Note that eo(sina) = e3(sina) = 0, e2(cosa) = e3(cosa) = 0, i.e., eola) = 
e3(a) = 0, or equivalently, 
(69) 
ex(o) = — cos aE (a) + sino.Ex(a) = 0, esla) = sina F(a) + cos aE3(a) = 0. 
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It follows that 


(70) F4(o) = 0, E3(a) = 0, 


which implies that the function o depends on only the variable v, i.e., a = a(v). 
Similarly, we can also conclude that the functions «1,0, fg and R = —qy, — q? 
depend on only v, i.e., kı = K1(v),0 = o (v), fo = falv) and R = —qy,—q? = R(v). 
In a similar way to those used computing the integrability data &1, fo and c in 
Eqs. (58)—(64), we can conclude the following 
(71) 

Kj = ki (v) = —sin? aqu, fo = falv) = — cos? aqu, e = o (v) = — sin a cos aqu, 


which means that the function q, depends on only the variable v, i.e., q, = v(v). 


Since e;(cos a) = e1(a3) = oa? = o sina, it is easy to see that 
2 D 


(72) a (v) 2 —o = sina cosaq,, 


where denote by a’ the first derivative of a along the geodesics curve v. Hereafter, 
we denote by f’, f" and f" the first, the second and the third derivatives of a 
function f along the geodesics curve v, respectively. Obviously, we have f' = 
ef) = y f” = eef) = 2d and f" = eiei(f) = TL Using Eq (72), one sees 


| .O0q | av) d(2a) 


(73) vv) m x - 


— 7A. du . E . 
Ov sinacosa  sin2a@ 


A further integration of the above gives 
(74) J ode + elu) = a = tena) + elu) 


which implies e?! = | tana(v)|?e??™. 

Consider the Coordinate transformation (u,v) > (u,v), u = f edu, DELE 
clearly, this is an isometric transformation M? — M? form M? to oneself. It 
follows that the Riemannian metric on M? can be also locally expressed as the 
form g = e??)qu? + dv? by reparametrization with e; = 2- and hence R = 
“Gy — q2 = —Q"(v) — Q? (v) on M?, where Q(v) = In|tano(v)|. It turns out 
that with respect to local coordinates (u, v, z), a product space M? x R can be also 
locally represented as (IR?, e?9 9 dy? + dv? + dz?) which is also obtained by using 
an isometric transformation ô : (RÌ, e4) qu? + dv? + dz?) > (R®, 22M dy? + 
dv? + dz?) with ó(u,v, z) = (u,v, z) = (f e?™du, v, z) from M? x R to oneself. 
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Combining these, (71),(72), and (73), we have 


kı = kv) = — sin? o(v)Q'(v) = — tan a o'(v), 
fo = fa(v) = — cos? aQ' = — cota a' (v), 
(75) o =o(v) = — sin a cos aQ’, 
a' (v) = —e = sina cos aQ'(v), 
Q(v) = In| tana(v)|, 
Q'(v) = ados: 


Finally, show that o satisfies the following ODE 
(76) a” sin a cos? a + cos o(sin? a + 3)oo" + sin a(2 cos? a + 3)o/? = 0. 


A straightforward computation using (75) gives 


E al? n H"H _3a'a! _ 2sinaa _ m 
(TT) = (ag + tana a"), Ki = Se ER tana a", 

100 ( a n 

2 = —(— z + cota a”). 


Noting that fı = k = 0, biharmonic equation (27) reduces to 
(78) Ak, — Ki{-K + fz} = 0, 
where KN = e;(fo) — f2. 
Using the results of Step 1, (20) and (36), one can easily compute the following 
Aki — kY( - KP + f?) 
(79) = €1€1(K1) — Veex(ra) — Vees(ri) — K1(—e1(f2) + 22) 
= K] — for, — kik — Kil fa + 2f2). 
We substitute (75), (77) and (79) into (78) and simplify the resulting equation 
to get (76). 
Step 2: With respect to new local coordinates, we show that the Riemannian 
submersion can be locally expressed as 
am: M? x R = (RÌ, 90d? + dv? + dz?) > (IR?, dv? + e? 0d9?), 
(u,v,z) = (v, 6) = (v, Ap Az + B), 
where the function Q(v) = In|tana(v)| and A = In|sina(v)| — In A depend on 
only v, and A > 0, B are constants. 


Let y : I —^ M? x R be an arbitrary integral curve of e1, which is geodesics 
since V4,e1 = 0. We parametrize y by y(v) = (u,v, z) by arc length parameter 
v. It follows from a known fact in |23, 24] that the curve 7 o y is a geodesic in 
N?, then we can parametrize the base space N? by rı = r1(¢,v) so that the v— 
curves are the integral curves of dz(ei) and the $— curves are the orthogonal 
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trajectories of the v— curves. Note that Gauss curvature of the base space KN = 
e(f2) — f2 = f5— f? = 307+ R cos? a by (23) and the 4th equation of (22), where 
R = —Q" — Q", so the Riemannian metric h on N? can be locally expressed as 
h = dw + e?™)d¢? with A = — fz, and hence —\"” — A? = KN = p — f? = 
3c? + Rcos? a. In addition, we use the 1st equation and the 4th equation of (75) 


1 . . . . 
to conclude X = $22 of which an integration gives 


(80) A =In|sina|+C; = ln | |+ Ci, 


Ce? 
where C, is a constant. 
From these, we see that the Riemannian submersion 7 : M? x R > (N?,h) can 
be locally expressed as the form m : M? x R = (R®,e?@™ dy? + dv? + dz?) > 
(R°, dv? + e? 9) q9?) with m(y, v, z) = (v, ¢), where ¢ is to be determined. 
We now determine ¢. Note that (e; = Ey,e2 = —cosaF, + sinaF3,e3; = 
sin aE, + cos BERT is an due frame adapted to m with e, being verti- 


cal and dz(ei) = 2. dr(e2) = e7 Es dr(e3) = 0. A further computation gives 
(81) 2 =dn(e:) = d«(2) = 24 #2, 
(82) F 
eg = = dr (e2) = dr(— cos a. E, + sinaE3) = — cos œe -9da(2- 4) t sin adn (2-) 
= (— cosae-9 Em + sinat) x, 
and 
(83) 0 = dr (e3) = dr(— sin aF; + cos a E3) = sin aae -?da(2- 5) + cos adn ( 2- =) 
= (sinae-9 Em + cos a 28) i. 
By comparing Ti sides of Eq. (81), one finds that 
do 
E =) 
Ov 


which means that the function o does not depend on v, i.e., ó = ó(p, z). 
Comparing coefficients of both sides of (82) and (83) separately, we get 


(84) — cos ae 9 — d + sin 420 =e, sin ae Bn + cos aoe eir. 
m Oz Oz 
We solve the above two equations to obtain 
(85) a = —e%)-®) cos a(v), ve =e sin a(v). 
Substituting the 5th equation of (75) and (80) into the above, we have 
Od o Od -8 
= = FA, Z = xe C = £A 
(86) T FA, 3,-— te l 
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where constant A > 0. 
Since ge = 0, the general solution of (86) can be computed as 


(87) d(u,2) = An d Az +B, 


where A > 0, B are constants. 
It is not difficult to check that if Q' = 0, i.e., œ is a constant, then the Riemannian 
submersion 7 is actually harmonic. 


Summarizing all results in the above cases we obtain the theorem. 


It is easy to note that a product space (IR?, e?» dx? + dy? - dz?) can be viewed 


as a twisted product space (IR? x2, R, dy? + dz? + e??™dx?). In particular, 


if p — p(y), the product space can be also considered as a warped product space. 
Therefore, Theorem 3.5 can be also stated as 


Corollary 3.6. Any proper biharonic Riemannian submersion v : R? x 2,4, R = 
(IR?, dy? + dz? + e?» dx?) — (N?,h) from a twisted (warped) product space is 
one of the following cases happens: 

(A): it exists only in R? X22.) R — R’. Moreover, it can be expressed as 


7 : (R2, dy? + dz? + e? 9 da?) > (R2, dy? + dz”), a(y,z,x) = (y, 2), 


and a function p(x,y) satisfies PDE (90) and p, 7 0, or, 
(B): it is a Riemannian submersion 


R?. dy? + dz? + e?» (dz?) > (R?, dy? + 8*9. dé?), 


(88) m ANLIRSRS 
my, z, v) = (y, $) = (y, FAT + Az + B), 
and a pde function p(x,y) = p(y) depends on only y and p(y) determined 


by a = arctan e") satisfying the ODE (35), where A 7 0, B are constants. 


Proof. The statement (A) is obvious by Theorem 3.5. For the statement (B), we 
first remind the readers to note that it is a convention to use coordinates (z, y, z) 
in place of coordinates (u,v, z) in the statement (b) of Theorem 3.5. Using the 
transformation u = x, v = y, z = z, a = arctan e" in the statement (b) of 
Theorem 3.5, we immediately have the statement (B). 


We can give many examples of proper biharmonic Riemannian submersions 


M? x R > R° by applying the following proposition stated as 
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Proposition 3.7. The Riemannian submersion 
(89) 7: M? x R= (R3, e?» dz? + dy? + dz?) — (R2, dy? + dz?) 
T(x, y, z) = (y, 2) 


is biharmonic if and only if 


(90) Pyyy + PyPyy + e P(e) (PrPzy = Peay) = 0. 

Furthermore, we have 

(a) the Riemannian submersion x is harmonic if and only if p = p(x), i.e., the 
function p depends on only x; 

(b) if p = ply), ie., the function p depends on only y, then Eq. (90) can be 
rewritten as the following equation 


(91) p" (y) + p'(y)p (y) = 0. 
Proof. Note that the frame (e1 — oe € = 2. e3 = e$) on (R°,g = 


cP) dx? + dy? + dz?) is adapted to the Riemannian submersion 7 with e3 being 
vertical, dz(e1) = €1 = » and dr(e2) = €2 = 2- form an orthonormal frame 
on the base space (R°, dy? + dz?). A straightforward computation gives the Lie 
brackets 


[e1, es] = —py es, [e2,e3] =0, [e1, e2] = 0, 
and hence the integrability data given by 
fi=fo=o =0, kı = —py, Ko = 0. 


m2 
For fı = fo = Ko = 0 and K R - 0, biharmonic Riemannian submersion equation 
(27) reduces to 
Ak, = 0. 


A further computation yields 
Pyyy + PyPyy + g "Pu (DzDzy Duy) = 0. 


From which, we obtain the proposition. 


Remark 3. (a): We would like to point out that the general solution of Eq. (91) 
can be obtained by theory of ODE. However, we are very interested in some 
particular solutions of Eq. (91). For example, we have the local solutions: 
p= C3y + Ca 
pln (Ci(sinh(Cay + C3))?) ; 
p= In (C1 (cos(Cay + C3))) : 
p = In (Cs(y + C3)?), 
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where C, 
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Cə > 0 and Cs, C4, ER. 


In addition to these we have the globally defined solution of Eq. (91) given by 


with C4, Co > 0, and C3 ER. 


= In (Ci(cosh(Coy + C3))?) , 


(b): The functions p = y(y)+ (2) give a family of solutions of Eq. (90), where 
the function y(y) is a solution of Eq. (90) and ó(x) depends on only z. 


(c): Since theory of the existence and uniqueness of solution of ordinary differ- 


ential equation, there is always a (local) solution o of Eq. (35) satisfying given 


initial condition. Therefore, for the solution o, we can conclude that there are 


many family of proper biharmonic Riemannian submersions defined by (34). It 


is easy to 


see that (34) can be rewritten as 


a: M? x R = (RP, tan? a(v)dy? + dv? + dz?) — (R?, dv? + 


n(p, v, z) = (v, $) = (v, FAW + Az + B), 


where A > 0, B are constant. 


By Remark 3, we immediately get the following 


Example 1. (i) (See also [34]) The Riemannian submersion 


a: H? x R= (R?, eda? + dy? + dz?) 


is a proper biharmonic map. 


(ii) The Riemannian submersion 


= 


a: H? xR = (R? xR, (cosh y)?da?--d?--dz?) > 


is a proper biharmonic map. 


R?, dy? + dz”), 


(iii) (See also [?, 6]) The Riemannian submersion 


mii 


R? x R, y*dz? + dy? + dz?) > (R2, dy? + dz”), n(x, 


is a proper biharmonic map, where y > 0. 


(iv) The Riemannian submersion 


T : (R? x 


R, (e*? cosh y)?da? + dy? + dz?) > 


( 


IR?, dy? + dz?), 


(R?, dy^--d2?), 


sin? alu sin a) d), 


nz, y, z) = 


y, z) = (y, 2) 


n(r,y,z) = 


is a proper biharmonic map, where the function ¢(a) depends on only z. 


m(x,y,z) = (y, 2) 


(y, z), 


(y, z) 
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